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On  Shock  Waves  in  Rock  Created  by 
Surface-  or  Near-to-Surface  Detonations 


by 

Leif  N.  Persen 


1 .Introduction. 

This  presentation  is  aimed  at  giving  a survey  of  the 
small  scale  experiments  performed  in  Norway  on  behalf  of  the  West- 
German  Defense  Ministry  to  start  a research  activity  in  the  field 
covered  by  the  title  of  this  report.  In  addition  the  data  from  the 
CENSE  experiments,  which  were  obtained  with  much  larger  charges, 
are  drawn  upon  to  see  if  conclusions  from  the  first  set  of  data 
could  be  substantiated.  The  small  scale  experiments  have  been  re- 
ported on  in  [l],  the  CENSE  experiment  is  covered  in  [2],  and  the 
data  from  the  latter  have  been  obtained  directly  from  the  record- 
ings. In  addition  results  from  some  earlier  experiments  are  drawn 
to  attention.  These  have  been  treated  separately  in  [3]. 

2 .General  remarks. 

The  general  philosophy  behind  the  small  scale  experi- 
ments referred  to  above  can  be  outlined  as  follows.  Usually  the 
rock  site  of  an  installation  has  been  tested  by  confined  explo- 
sions. In  this  way  the  properties  of  the  rock  site  as  a transmitt- 
ing medium  may  be  considered  known,  for  instance  in  the  way  descri- 
bed in  [3].  Consequently  the  creation  of  a shock  wave  and  its  trans- 
mission through  the  rock  may  be  said  to  be  well  known,  provided  the 
magnitude  of  the  creating  charge  is  known,  and  that  it  is  detonated 
fully  confined. 

In  the  case  of  a shelter  of  some  kind,  the  charge 
creating  the  shock  wave  is  usually  considered  to  be  a surface 
charge  of  some  sort.  Detonation  is  supposed  to  take  place  above 
ground  level,  on  the  ground  or  at  a certain  penetation  depth  be- 
neath the  surface.  The  problem  then  arises  how  to  make  predictions 
for  the  propagation  of  the  created  shock  wave  based  only  on  the 
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knowledge  of  the  fully  confined  case.  One  way  of  answering  this 
question  is  to  perform  exeriments  whereby  the  position  of  the 
charge  above  (or  beneath)  the  surface  is  changed  from  one  detona- 
tion to  another,  each  time  observing  the  created  shock  wave.  This 
idea  was  more  or  less  the  basis  for  the  CENSE  experiments,  and 
the  result  came  out  as  diagrams  of  uie  "containment  factors"  for 
different  quantities. 

The  small  scale  experiments  were  based  on  a somewhat 
more  sophisticated  philosophy.  As  shown  in  [3],  pp.192-19^,  a possi- 
bility exists  to  draw  direct  conclusions  for  socalled  "half  burried" 
cases  from  knowledge  of  the  fully  contained  case.  It  was  therefore 
thought  sufficient  to  relate  cases  with  varying  heights  of  the 
detonating  charge  above  the  surface  to  the  "half  burried"  case.  In 
this  way  the  problem  was  split  in  two  "stages",  each  stage  being 
basically  different  from  the  other  as  the  physics  indicate.  For 
positions  of  the  charge  on  or  close  to  but  beneath  the  surface 
cratering  becomes  a problem.  For  positions  of  the  charge  above  the 
surface  cratering  does  not  necessarily  take  place.  It  seems  realis- 
tic to  treat  these  two  "stages"  separately. 

3. The  experiment  Vikersund L 

A small  scale  experiment  were  performed  in  VIKERSUND, 
Norway.  Two  different  geometries  were  used.  The  first  one  is  shown 
in  Fig.l  . A natural  precipice  with  a rather  un-eroded  surface  was 
located.  In  front  of  or  at  the  vertical  wall  the  charges  were  de- 
tonated at  A.  From  the  top  of  the  precipice  the  holes  1-9  were 
drilled  so  that  the  pick-ups  used  for  measuring  strain  or  particle 
acceleration  as  the  shock  wave  was  transmitted  could  be  positioned 
in  the  rock.  All  measuring  devices  were  placed  on  the  perpendicular 
to  the  surface  at  A.  A separate  hole  was  drilled  so  that  a confined 
charge  could  be  detonated  in  B.  This  made  it  possible  to  monitor 
the  shock  wave  as  it  travelled  in  both  directions  between  A and  B, 
and  thereby  eventual  systematical  deviations  could  be  detected. 

This  part  of  the  experiment  was  more  or  less  a repetition  of  the 
experiment  described  in  [3],  pp  175-19^. 
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Fig.l.  Geometry  of  the  "row" -experiment 
at  VIKERSUND , Norway. 

In  the  second  geometry  used  at  VIKERSUND,  the  tech- 
nique of  fixing  strain  gauge  directly  to  the  rock  was  used.  A hole 
was  drilled  normal  to  the  surface  of  the  rock  as  shown  in  Fig. 2 . 
Strain  gauges  were  fixed  to  the  rock  surface  at  different  distances 
in  the  hole  from  the  surface.  The  necessary  wiring  was  led  away  in 
a sloping  hole  to  prevent  it  from  being  broken  by  the  explosion. 

The  hole  is  then  filled  with  concrete,  and  it  was  hoped  that  the 
signals  would  give  a direct  measurement  of  the  strain  in  the  rock 
during  the  passage  of  the  shock  wave.  The  shock  wave  was  created 


Fig. 2.  Geometry  of  the  "hole"-experiment  at  VIKERSUND y Norway. 

by  letting  a charge  be  detonated  as  shown,  whereby  the  distance  h 
and  the  charge  magnitude  W were  varied  from  case  to  case.  The  re- 
sults from  these  experiments  will  be  treated  separately. 
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4. The  "row"-experiment,  VIKERSUND. 

In  the  "row"-experiment  only  charges  in  contact  with 
the  surface  of  the  rock  at  A in  Fig. I were  detonated.  Consequently 
an  increasing  crater  was  developing  at  A as  the  experiment  progress- 
ed. Because  the  magnitude  W of  the  charge  was  increased  gradually, 
one  operated  almost  automatically  with  "half  hurried"  charges. 

The  shots  were  numbered  consequtively,  and  Table  I 
shows  how  the  charge  magnitude  was  changed.  The  table  also  give 
the  maximum  amplitude  of  the  shock  wave  measured  in  [y-strain]  at 
each  of  the  nine  pickups.  These  details  are  given  to  illustrate 
the  repeatability  of  the  measurements.  As  shown  the  repeatability 
is  not  as  good  as  one  could  wish,  and  each  signal  from  each  pickup 
was  therefore  scrutinized  for  "misbeaviour"  which  might  indicate 
decreased  reliability.  It  was  found  that  only  the  shots  A4,A5,A6, 
A11,A12  and  A14  could  be  used  with  some  confidence.  The  fully  con- 


TABLE I 

Maximum  amplitude  A^ [y-s train] 
measured  at  pickup  No: 

W 

[kg] 

t-1 

2 

3 

4 

5 

6 

7 

8 

9 

B 

1 

0,5 

61 

83 

131 

157 

183 

256 

610 

792 

B 

2.0 

58 

71 

71 

- 

- 

214 

490 

875 

A 

0.5 

792 

238 

- 

- 

- 

_ 

- 

- 

A 

0,5 

582 

155 

22 

21 

11 

9 

— 

- 

- 

A 

0,5 

629 

190 

36 

34 

15 

12 

- 

22 

17 

A 

2,0 

1780 

547 

155 

131 

45 

42 

21 

(32) 

29 

A 

7 

2,0 

1466 

309 

(27) 

42 

20 

16 

- 

18 

12 

A 

8 

2,0 

- 

285 

95 

98 

46 

36 

- 

41 

24 

B 

9 

2,0 

23 

36 

60 

131 

78 

131 

tm 

626 

933 

A 

10 

5,0 

1982 

1368 

268 

- 

- 

- 

- 

128 

63 

A 

11 

5,0 

_ 

893 

190 

238 

76 

58 

59 

85 

51 

A 

12 

15,0 

- 

2320 

690 

547 

218 

204 

140 

321 

168 

A 

13 

15,0 

- 

2507 

489 

319 

171 

148 

154 

294 

139 

A 

14 

30,0 

- 

4160 

1082 

786 

278 

296 

357 

514 

274 

A 

42 

30,0 

— 

4047 

1071 

785 

171 

185 

330 

350 

245 

tained  shot  Bl,  B2  and  B9  will  be  treated  senarately,  as  these  are 
giving  the  reference  frame  for  the  wave  propagation  for  this  par- 


ticular  rock,  (the  norwegian  Gneis  ). 

The  arrival  times  for  each  signal  at  each  pickup 
was  used  to  determine  the  signal  velocity  in  the  rock.  Because  the 
positions  of  the  pickups  are  known,  the  arrival  times  t ^ may  be 


Fig.  3.  Arrival  times  'plotted  as  funotion  of  the  distance 

plotted  as  a function  of  the  distance  d^  of  each  pickups  from  the 
point  of  detonation  A,  as  shown  in  Fig.  3.  The  slope  of  the  straight 

line  through  the  points  will  then  give  the  signal  velocity  c of 

s 

the  rock.  Now  there  are  two  possibilities  for  defining  this  line: 


either:  d - d + o t 

o s 


d - a t 
s 


Linear  regression  gives: 


either:  o ~ 6327  [m/s]  , d = 0.27  [m] 
s o 

or:  c — 6030  [m/s] 
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The  two  cases  are  illustrated  in  Fig. 3 by  a dashed  line  and  a 
fully  drawn  line  respectively.  The  difference  between  the  to 

results  amounts  to  5%,  and  the  absolute  value  of  a is  rather 

s 

high.  This  rather  trivial  determination  of  the  signal  velocity 
has  been  given  in  detail  here  because  it  indicates  the  quality 
of  the  measurements. 

The  data  for  the  "row"-experiment  conslv  ? of  the 
maximum  amplitude  /l.  measured  at  pickup  No./,  the  dispersi  n C\ 
measured  as  the  maximum  slope  of  the  signal  in  it  first  pi  -ise 
at  each  pickup,  and  the  distance  from  the  explosion  to  the 

pickup.  In  addition  the  calibration  constant  k.  of  each  pickup 
may  be  considered  as  input  data.  These  are  now  treated  as  shown 
in  [3]  to  obtain  the  attenuation  curves  for  the  maximum  amplitude 
and  the  dispersion.  The  whole  procedure  is  given  in  Tables  II  and 
III,  where  the  determination  of  the  shot  factors  3,  the  corrected 

values  of  the  non-dimensional  distances  L.  as  well  as  of  the 

aorv 

charactreistic  lenghts  a for  each  shot  are  also  shown.  The  disper- 
sion data  have  been  made  dimensionless  using 

c - 6250  [m/s]  , p = A*  (4.5) 

s ^00 

The  attenuation  curves  will  assymptotically  be  ex- 
pressible as 

A = A*£^  (for  the  max.  ampl.  ) 
o 

A 

C = C*6,  a (for  the  dispersion) 

The  evaluation  delivers  the  values  of  A*.C*.  A and  A , with  the 

o o a 

added  information  of  the  standard  deviation  m given  in  percent  of 
the  measured  values. 

As  mentioned  some  of  the  data  obtained  were  rejected 
leaving  only  data  from  6 of  the  shots  to  be  considered.  To  make 
sure  that  one  has  not  introduced  any  irrelevant  regularity  by  this 
selection,  the  shots  A12  and  A14  with  the  larger  charges  were  also 


(4.6) 

(4.7) 


18.11  8.26  ; 3.C1S06 

21.29  6.90  3.15001 


TABLE  IV  TABLE  v 


> ao  30  40  50  100  10  to  30  40  60 

« < 

Fig.  4.  The  resulting  attenuation  curves  for  the  maximum  amplitude  A 
and  for  the  non-dimensional  "dispersion"  C .a/p  c k. 

t O S Is 

(Half  hurried  charges) 
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neglected  and  the  evaluation  repeated  as  shown  in  the  Tables  IV 
and  V.  Inspection  of  the  data  reveals  that  only  data  from  the 
pickups  No. 2, 3,4, 5 and  6 can  be  relied  upon.  Those  data  are  marked 
with  a star  (*)  in  the  tables. 

The  results  for  the  six  shots  as  well  as  for  the  four 
shots  are  as  follows: 


For  the  6 shots: 
(Tables  II  and  III) 


For  the  4 shots: 
(Tables  IV  and  V) 


A*  = 41576  [atm] 
o 

A - -2.15234 
m = ± 17.8 % 


A*  = 43765  [atm] 
o 

A = -2.16719 
m = ±18.5% 


(4.8) 


C = 10.88044 
o 

A = -3.11369 
a 

m = ±24.  7% 


C = 11.83964 
o 

A = -3. 14128 
a 

m = ±23. 6% 


(4.9) 


For  all  practical  purposes  these  results  may  be  considered  as 
identical . 

Tables  IV  and  V show  the  legend  used  when  the  result  is 
plotted  in  Fig. 4 . 


Unfortunately  the  data  from  the  fully  confined  charges 
detonated  in  B in  Fig.l  are  not  reliable.  The  reason  for  rejecting 
data  from  the  pickups  No. 7, 8 and  9 in  the  evaluation  of  the  sur- 
face charges  was,  that  the  signals  showed  a behaviour  which  could 
be  interpreted  as  caused  by  a major  crack  in  the  rock  between  pick- 
ups No. 6 and  7.  In  this  case  the  shock  wave  travels  in  the  oppo- 
site direction,  and  one  would  expect  only  data  from  the  pickups 
No.  9,8,  and  7 to  be  relevant.  The  data  are  plotted  in  Fig. 5 and 
it  is  observed  that  the  wave  from  the  0.5  kg  charge  is  damped  less 
than  the  wave  from  the  2.0  kg  charge.  Consequently  the  smaller 
charge  creates  a larger  maximum  amplitude  in  the  shock  wave 
after  a certain  distance  than  the  larger  one,  a result  which  must 
be  rejected  on  physical  grounds. 
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Fig.  5.  Data  from  shots  Bl,  B2  and  B9  with  the 

attenuation  curves  derived  from  the  data 
of  Bl  and  B9  seperately . 


Difficulties  occurred  with  the  charge  in  shot  B2.  The  data 
from  this  shot  are  therefore  neglected.  The  dashed  lines  in  Fig.  5 
show  the  attenuation  curves  given  by  the  two  shots  Bl  and  B9  seperately. 
The  deviation  is  so  great  that  further  attempts  at  drawing  meaningful 
conclusions  from  these  data  are  thought  futile.  This  means  that  the 
previously  established  correlation  between  the  effect  of  a half  hurried 
charge  as  compared  with  a fully  contained  charge  could  not  be  checked. 


5.  The  “hole"-experimentt  VIKERSUND. 

The  geometry  and  the  purpose  of  the  "hole"-experiment  at 
Vikersund  has  already  been  treated  and  shown  in  Fig. 2 . Fig. 6 shows 
the  17  strain  gauges  and  their  positions  relative  to  the  surface  of 


Fig. 6 Positions  of  the  strain  gauges  in  the  "hole "-experiment 
at  VIKERSUND  given  in  centimeters  from  the  surface. 

the  rock.  It  is  noticed  that  the  gauges  5 and  6,  11  and  12,  and 
14  and  15  are  placed  opposite  each  other  at  the  same  distance  from 

TABLE  VI 
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the  surface.  The  maximum  amplitudes  of  the  shock  waves  measured  in 
[y-strain]  at  each  strain  gauge  as  well  as  the  specifications  of 
each  shot  (distance  above  the  surface  h and  charge  magnitude  U ) 
are  given  in  TABLE  VI. 

It  is  observed  that  the  instrumentation  of  this  experi- 
ment is  different  from  the  one  used  in  the  "row"-experiment . The 
question  of  whether  or  not  the  type  of  instrumentation  influences 
the  data  such  that  conclusions  drawn  on  the  physical  behaviour  really 
reflects  instrument-induced  errors  becomes  important.  It  is  there- 
fore imperative  somehow  to  connect  the  results  from  this  experi- 
ment with  that  of  the  "row"-experiment . This  occurs  first  through 
the  determination  of  the  signal  velocity  of  the  rock.  In  Fig. 7 are 
the  arrival  times  t of  the  signals  at  the  different  pickups  and 
the  distance  d travelled  plotted  against  each  other  for  the  shots 
16,  18,  20,37  and  38.  The  data  exhibit  a comparatively  large  scatter. 


O lOO  200  300 

‘aO'Vl  - 

Fig. 7 . Arrival  time  t of  the  shook  wave  at  the  pickups  related 
to  the  distance  d travelled.  "Hole"-experiment,VIKERSUND: 
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but  the  straight  line  through  these  points,  determined  by  linear 

regression,  gives  through  its  slope  the  signal  velocity  a as 

s 

follows : 


c - 6210  [m/s] 
3 


b - 5.18  [m] 


(5.1) 


where  b is  the  distance  at  which  t = 0,  leaving  an  impression  of 
the  accuracy  of  the  measurements.  The  agreement  between  the  results 
(4.3),  (4.4)  and  (5.1)  is  satisfactory. 

The  only  data  in  Table  VI  which  may  be  compared  to  the 
results  obtained  in  the  "row"-experiment  are  the  ones  from  the 
shots  40  and  41.  Only  in  these  two  cases  did  one  have  the  charges 
in  contact  with  the  rock's  surface.  The  data  are  repeated  in 
Table  VII  which  also  gives  the  legend  for  the  points  in  Fig. 8 . 

The  straight  line  through  the  data 
TABLE  VII  points  will  have  a slope  A which 

can  be  determined  by  linear  re- 
gression: 


A - -1.9779  (5.2) 

This  value  should  be  compared  with 
the  values  obtained  in  (4.8).  A 
better  idea  of  the  agreement  be- 
tween the  two  cases  is  however 
obtained  by  comparing  the  results 
in  an  attenuation  diagram.  Fig. 8 
shows  this  where  the  fully  drawn 
line  represents  the  attenuation 
curve  determined  from  the  data 
shown,  and  the  dashed  line  repre- 
sents the  attenuation  curve  ob- 
tained from  the  results  of  the 
"row"-experiment  with  a = 0.06635  [m] 
for  a charge  W = 0.5  [kg]  and  with  a mean  value  of  k = 3. 1 [p-str/atm] 
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Fig.  8 . Maximum  amplitudes  A.  ( "hole”-experiment) 
by  the  shots  40  and  \l. 


for  the  calibration  constant.  Considering  the  fact  that  one  is 
in  this  case  comparing  results  from  experiments  with  two  entirely 
different  types  of  pickups,  the  agreement  must  be  considered  very 
satisfactory. 

So  far  the  experiments  have  dealt  with  situations  for 
which  at  least  some  sort  of  theoretical  approach  could  be  en- 
visaged. For  cases  with  charges  detonated  in  the  air  above  the 
surface,  no  theoretical  method  known  to  the  author  exists  by 
means  of  which  one  can  relate  the  shock  wave  induced  in  the  ground 
to  that  which  the  same  charge  would  have  caused  when  detonated 
completely  confined.  The  idea  behind  the  experiments  was  to  do 
this  experimentally. 
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The  free  hanging  charges  will  induce  shock  waves  in  the 
rock  which  will  be  influenced  among  other  things  by  the  distance  h 
above  the  surface  at  which  the  charge  is  detonated.  This  distance 
may  be  made  dimensionless  by  scaling  it  against  the  charge  mag- 
nitude W , i.e.  h*  = h/\f\T  . The  experiment  was  planned  such  that 
the  non-dimensional  distance  h+  was  repeated  for  each  charge  mag- 
nitude. One  observes  that  this  means  equivalence  in  the  way  in 
which  distances  are  scaled  in  the  rock  as  well  as  above  its  sur- 
face . 

Before  examining  the  experimental  data  one  may  contemp- 
late what  kind  of  relationship  one  may  anticipate.  If  one  uses 
the  attenuation  curve  for  the  maximum  amplitude  of  the  shock  wave 
in  the  case  of  fully  contained  explosions  as  a guide,  one  will 
expect  an  attenuation  curve  as  shown  in  Fig. 9 , where  the  maxi- 


Fig.9  . Attenuation  curve  for  the  maximum  amplitude  A 


mum  amplitude  A is  plotted  as  a function  of  the  distance  d travel- 
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led  by  the  shock  wave.  The  diagram  is  in  doubly  logarithmic  scale 
and  consequently  the  curve  ought  to  exhibit  an  assymptotic  be- 
haviour as  a straight  line  as  shown  by  the  dashed  line. 

One  may  be  more  specific  at  this  point  . One  cannot  off- 
hand know  how  the  attenuation  curve  is  influenced  by  the  rock 
type,  the  type  of  explosive  used,  the  dimensionless  duration  of 
the  input  pulse  and/or  the  elevation  h above  the  surface.  One  can 
however  assume  as  a working  hypothesis  that 

a)  the  distances  may  be  scaled  with  VvT  , 

b)  the  dimensionless  duration  x is  the  same  in  all 
cases}  because  among  other  tRings  the  same  explo- 
sive is  used  , 

c)  the  data,  which  here  consists  of  measured  maximum 
amplitudes  of  the  shock  wave,  will  depend  only  on 
one  parameter,  say  h*. 

Such  a hypothesis  corresponds  to  the  expected  behaviour  of  the 
attenuation  curve  for  a contained  explosion. 

This  hypothesis  is  now  examined  in  the  two  Tables  VIII 
and  IX.  First  Table  VIII  is  arranged  in  such  a way  that  the  three 

cases  of  constant  h + appear  separately.  One  then  computes  the  ’>] 

ratios  between  the  measured  maximum  amplitudes  at  each  pickup. 

0 5 y ^20  and  A50  are  the  maximum  amplitudes  measured  with 
W =0.5[kg]  , W = 2.0[kg]  and  W = 5.0[kg]  respectively.]  It  is 
observed  that  the  ratios  remain  fairly  constant  for  larger  distan- 
ces. The  average  value  of  the  ratios  neglecting  those  closest 
to  the  charge  and  appearing  in  paranthesis  ()  in  the  table  is 
given. 

Secondly  Table  IX  is  arranged  such  that  the  three  cases 
of  different  charge  magnitudes  appear  separately.  Again  the  ratios 
are  computed  and  again  the  average  values  of  the  ratios  are  given 
for  the  pickups  farthest  away  from  the  surface.  This  exercise 
indicates  that  the  data  should  be  subjected  to  a scaling  procedure 
which  is  exhibited  in  Fig. 10  . On  the  transparent  page  1 the 
data  from  the  0.5  kg- charges  at  different  distances  h are  gathered. 

The  scaling  occurs  vertically  and  is  proportional  to  1/h On 
page  2 the  same  has  been  done  for  the  2.0  k^-charges,  but  here  a 
horizontal  scaling  occurs  in  accordance  with  point  a)  above.  On 
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page  3 the  same  procedure  is  repeated  for  the  data  from  the  6.0  kg- 
charges.  Because  the  pages  are  transparent  the  figure  reveals  how 
nicely  the  data  seem  to  follow  the  scaling  laws  outlined  above. 

One  should  at  this  point  add  the  following  remarks: 

The  ratios  computed  may  to  some  extent  be  independent  of  eventual 
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systematic  errors  in  the  data.  An  example  is  the  data  from  pick- 
up no.  14  which  evidently  are  erroneous.  Still  the  ratio  computed 
from  them  is  not  unrealistic.  This  means  further  that  even  though 
the  data  give  ratios  with  fair  degree  of  accuracy  ( <10% ) they 
may  still  exhibit  a rather  large  scatter.  Fig. 10  shows  how  this 
is  the  case  for  data  obtained  close  to  the  surface. 
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It  must  be  remarked  that  the  data  for  the  case  of 
W - 0.5  kg  and  h - 20  am  do  not  seem  to  fit  in  with  the  proposed 
scaling  procedure.  It  is  therefore  plotted  in  with  its  base  line 
dashed . 

Finally  it  should  be  remarked  that  the  exercise  so  far 
has  only  given  a scaling  procedure  for  cases  of  free  hanging  charges. 
If  it  is  extended  to  the  case  of  the  charge  in  contact  with  the 
surface,  ( h + - 1 ),  one  would  have  been  able  to  connect  the  pre- 
sent result  to  the  case  of  half  burried  charges  and  thus, through 
the  results  of  [3],  also  to  the  case  of  a confined  charge.  In 
Fig. 10  the  data  from  shots  40  and  41  are  plotted  on  page  4 realizing 
that  these  data  may  be  regarded  as  a case  for  which  h * - 2, where- 
as the  other  cases  give  values  of  this  parameter  as  shown  in 
Table  VIII  . The  fact  that  these  data  fit  in  so  nicely  with  the 
rest  may  be  taken  as  a very  strong  indication  that  the  scaling 
procedure  outlined  here  may  be  acceptable. 


6.  Conclusions  from  the  VIKERSUND  experiments. 


It  is  not  easy  to  draw  general  conclusions  from  experi- 
ments within  a limited  range.  Keeping  this  in  mind  one  is  however 
encouraged  by  the  apparent  consistency  in  the  data  and  an  attempt 
may  be  made  to  generalize  the  results  as  follows: 


When  comparing  two  cases  I and  II  where  the  charges 
Wj  and  W are  being  detonated  in  the  same  distance 
h*  above  the  surface  of  the  rock,  one  will  find  that 
the  same  maximum  amplitude  in  the  shock  wave  is  created 
at  distances  dj  and  d respectively,  whereby 


(6.1) 


(Maximum  amplitude  measured  in  [y-strain]) 


This  is  to  be  considered  the  first  part  of  a "scaling  law"  or  a 
"model  law".  The  non-dimensional  distance  h*  used  here  is  defined 


h * = h/R 


R = Vsw/4-nQ 


(6.2) 


where  h is  the  elevation  above  the  surface  where  the  detonation 
takes  place,  and  R is  the  radius  of  the  equivalent  sperical  charge 
with  p as  the  density  of  the  explosive.. 

The  second  part  of  the  "scaling  law"  may  be  formulated 
as  follows: 


When  comparing  two  cases  I and  II  where  the  same  charge 

W is  detonated  at  the  non-dimensional  elevations  T and 

J I 

h*  respectively,  one  will  find  that  in  the  same  non- 
dimensional  distance  from  the  surface  the  shock  wave 
will  exhibit  maximum  amplitudes  S ^ and  respectively 
which  will  be  related  to  each  other  through 


II 


*,II 


K 


(6.3) 


If  this  result  can  be  substantiated  and  confirmed  by  further 
experimental  evidence,  one  has  succeeded  in  relating  at  least  one 
important  case  of  air  blast  induced  shock  waves  in  rock  to  the 
case  of  a confined  explosion.  For  cases  where  crater  building 
is  a major  factor  in  the  process,  reference  is  being  made  to 
the  handling  of  this  problem  in  [3]  • 
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7.  The  CENSE  experiments. 

In  connection  with  the  VIKERSUND  experiments, 
which  were  conducted  on  a small  scale,  it  may  be  of  interest  to 
examine  the  results  of  the  CENSE  experiments.  These  were  conducted 
with  charges  which  were  orders  of  magnitude  greater.  The  results 
from  these  experiments  are  discussed  in  [2]  and  the  lay-out  of 
the  experiments  is  shown  in  Fig. 11  . The  numbers  indicate  the  7 
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Fig. 11.  The  experimental  lay-out  of  the  CENSE  experiments. 

different  events,  where  for  each  event  the  position  of  the  deton- 
ating charge  is  changed.  In  this  way  one  covers  the  whole  range  of 
positions  above,  on,  near  to  and  below  the  surface.  For  each  event 
recordings  of  the  induced  shock  wave  in  the  rock  were  made 
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on  a vertical  ray  through  the  center  of  the  charge  as  indicated 
through  the  vertical  instrument  array.  In  case  7 also  a horizon- 
tal instrument  array  was  in  operation.  Other  measurements  were 
also  made,  but  these  are  so  far  of  marginal  interest  in  the  pre- 
sent context. 

It  should  be  mentioned  that  the  results  of 
these  experiments,  as  they  appear  in  [2],  are  given  in  terms  of 
coupling  factors.  As  an  example: the  ground  shock  factor  F is 
defined  as 


F 

u 


oeak  horizontal  particle  velocity 


peak  radial  velocity  for  full  containment  (DoB  = 7R J 


It  is  stressed  that  the  ground  shock  factor  is  determine  ’ from 


measurements  along  the  near  surface  instrument  radial,  and  that 
data  from  the  measurements  directly  beneath  the  explosive  are 


not  included  in  the  analysis.  This  is  in  contrast  to  the  present 
examination,  where  only  the  data  from  the  vertical  instrument 
array  will  be  used.  (In  event  7 also  the  horizontal  array  will 


be  considered.) 


8.  Peak  velocity  data  from  the  CENSE  I experiment. 


The  data  from  CENSE  I are  presented  in  the 
form  of  recorded  particle  velocity  as  function  of  time  as 
examplified  in  Fig. 12  . The  peak  particle  velocity  is  rather 


Fig.  12.  Recorded  particle  velocity  at  a distance  of  20  feet 
in  Event  S. 
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easily  picked  out  from  such  recordings,  and  this  quantity  is 
used  to  characterize  the  shock  waves  as  done  also  in  the  first 
report  [2].  These  data  are  listed  in  Table  X as  the  maximum 

amplitude  A.  measured  at  pickup  no.  i in  a distance  d.  from  the 

x *• 


explosion. 


Fig.  12.  Data  from  event  7. 


One  may  start  examining 
these  data  by  looking  at  the  fully 
contained  case,  event  7.  This  is 
to  be  a reference  case  and  it  is 
imperative  that  this  case  is  well 
established.  If  the  data  are  plot- 
ted as  done  in  Fig. 12  , it  becomes 
however  apparent  that  the  attenua- 
tion of  the  maximum  amplitude  of 
the.  shock  wave  is  much  greater  in 
the  vertical  than  in  the  horizon- 
tal direction.  This  may  indicate 
that  the  rock  (sandstone)  exhibits 
different  properties  as  a wave- 
transmitting  medium  in  the  two  di- 
rections, i.e.  the  sandstone  is 
"layered".  Or  it  may  reflect  the 
fact  that  one  or  more  of  the  pick- 
ups used  exhibit  systematic  errors. 
Because  repeated  experiments  were 
not  undertaken,  and  no  other  infor- 
mation is  available,  the  question 
of  which  of  the  two  possibilities 
is  really  the  case  must  remain  un- 
resolved . 

Both  data  from  the  horizon 
tal  and  from  the  vertical  array  in 
event  7 record  what  may  be  assumed 
to  be  a sherical  wave.  The  idea  be- 
hind the  handling  of  the  data  from 
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the  charges  placed  on  or  near  to  the  surface  is  to  assume  that 
the  situation  on  the  vertical  through  the  charge  can  be  compared 
with  the  confined  case.  Now,  the  confined  case  appears  to  give 
ambiguous  informations,  and  the  reference  case  thus  becomes  un- 
certain . 

In  [2]  the  ground  shock  factor  F is  defined 

u 

in  accordance  with  (7.1)  where  the  data  from  the  horizontal  ar- 
ray are  used  for  the  reference  case.  If  a straight  line  is  fit- 
ted through  the  data  points  (a)  from  the  horizontal  array  in 
Fig. 12  , an  approximation  to  the  attenuation  curve  is  obtained 
which  may  be  expressed  as 

A = 10752  d.2'03052  [ft/sec ) (8.1) 

% 

where  A is  the  peak  particle  velocity  of  the  shock  wave.  The 
mean  deviation  of  the  data  points  from  this  line  is  ±35-1%. 

The  slope  A of  the  attenuation  curve  determined 
from  the  data  of  the  horizontal  array  above  is  equal  to  the  ex- 
ponent -2.03053  in  the  expression  (8.1).  The  corresponding  slope 
determined  from  the  data  of  the  vertical  array  (a)  in  Fig. 12 
is  much  greater.  Even  if  the  point  originating  from  pickup  no. 5 
is  neglected  as  a stray  point,  the  value  of  the  slope  would 
be 

A - -2.66817  (8.2) 

and  the  mean  deviation  of  the  points  from  the  straight  line 
shown  in  Fig. 12  would  be  ±32.6%.  The  difference  in  magnitude 
of  the  slopes  as  well  as  the  vertical  position  of  the  straight 
lines  representing  the  attenuation  curves  in  Fig. 12  is  clearly 
brought  out  by  the  plot.  One  has  however  a possibility  to  check 
this  result  by  drawing  Fig. 3 of  the  report  [2]  to  attention. 

This  figure  is  reproduced  in  Fig. 13  with  the  original  figure 
caption.  It  is  seen  that  the  different  cases  give  lines  of  ap- 
proximately the  same  slope  A -1.55*  This  indicates  that  the 
present  investigation  may  be  based  on  data  which  differ  from 


28  - 


Fig. IS.  Replot  of  Fig.  S in  reference  [2]  with 
the  caption: "Peak  particle  velocity 
versus  scaled  range  as  a function  of 
charge  containment". 

those  used  in  reference  [2] . It  should  however  be  noted,  that  the 
magnitude  of  the  slope  in  Fig. 13  corresponds  well  to  the  slope 
A = -1 .65  suggested  by  Fred  M.  Sauer  in  Nuclear  Geoplosics  [2]  for 
the  socalled  composite  attenuation  curve,  and  that  the  slopes 
obtained  from  the  present  data  in  (8.1)  and  (8.2)  seem  unusually 
high  by  comparison. 

In  spite  of  the  discrepancy  which  seem  to  have 
been  discovered,  it  may  well  be  that  the  data  give  adequate  re- 
sults for  the  ground  shock  factor  F as  defined  in  (7.1).  This 

u 

factor  is  therefore  computed  for  each  shot  and  each  distance  in 

Table  X.  Also  the  factor  F*  with  the  data  from  event  5 as  refer- 

u 

ence  has  been  calculated  in  the  same  table.  Both  cases  reveal 
that  no  obvious  trend  can  be  found  in  these  results. 


TABLE  XI 


Event 

3 

1 

■ 

A . 

X 

[ ft/ sec ] 

B. 

3 

B A . 

3 1 

[ ft/sec] 

1 

3 

40 

.463 

1.82616 

.846 

4 

50 

.219 

. Uoc 

O 

5 

60 

.147 

.268 

2 

1 

20 

1.767 

2.22640 

2 

30 

• 514 

3 

40 

( .432) 

.962 

4 

50 

.170 

.376 

• 

5 

60 

. 112 

.249 

3 

1 

20 

9.32 

1.53131 

14.27 

2 

30 

3. 1 84 

4.876 

3 

4o 

.844 

1.292 

4 

50 

. 192 

.294 

□ 

5 

6o 

.156 

.239 

4 

1 

20 

12.85 

.95670 

12.29 

2 

30 

6.077 

5.81 

3 

4o 

.836 

.800 

4 

50 

.620 

.593 

■ 

5 

60 

.200 

. 191 

5 

1 

20 

22.64 

.58801 

13.313 

2 

30 

5.10 

2.999 

3 

Lo 

1.29 

• 759 

4 

50 

.828 

.487 

O 

5 

60 

.418 

.246 

6 

1 

20 

131.00 

.28552 

37.40 

2 

30 

15.50 

4.425 

3 

40 

2.203 

.629 

4 

50 

2.452 

.700 

♦ 

5 

60 

.722 

.206 

In  spite  of  this  rather  discouraging  result,  one  may 
attempt  to  establish  the  ground  shock  factors  by  introduction  of  the 

"shot  factors"  |3  . as  shown  in  [3]  and  then  determine  these  by  means 

3 

of  the  method  of  the  least  squares.  The  result  is  shown  in  Table  XI 
where  the  data  from  Event  7 have  been  disregarded  entirely.  The  result 

ought  to  be  such  that  the  data  (d.  , (3  .A  .) , when  exhibited  in  an 

i 0 ** 

attenuation  diagram  , would  gather  around  a straight  line  with  the 
slope  A.  In  the  present  case  the  calculation  is  based  on  the  data 
from  the  pickups  3 ,4  and  5 only,  and  the  result  is  shown  in  Fig. 14  . 
The  slope  A of  the  straight  line  and  the  mean  deviation  m of  the 
points  from  this  line  will  be 
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Fig.  14  . Corrected,  data  3 -A.  plotted  as  function  of 
the  distance  a . . ^ ^The  straight  line  repre- 
sents the  attenuation  curve  for  particle 
ve loci ty . 


X = -3.21792  , m-  ±27. 07c  (8.3) 

It  should  be  noted  that  the  scatter  indicated  by  m only  reflects 
the  deviation  of  the  data  from  the  pickups  3 , 4 and  5.  Fig. 14 
reveals  a much  greater  scatter  as  far  as  the  data  from  pickups 
no.  1 and  2 are  concerned. 


M 
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Fig. 15  . Corrected  data  g ,B . plotted  as  fun  tion  of 
the  distance  d . 'I  1 The  straight  line  repr 
sents  the  attenuation  curve  for  particle 
ve  loci  ty . 


It  must  be  mentioned  that  the  data  can  be  obtained 
from  a second  set  of  curves,  which  is  given  with  the  first  one. 

This  set  is  not  as  complete  as  the  first  one,  the  peak  particle 
velocities  obtained  from  it  are  denoted  by  , and  they  are 
given  in  Table  XII,  where  for  comparison  the  data  A.  are  repeated. 
Missing  information  is  denoted  by  a star  (*)  in  the  table,  and 
illegible  signals  with  a bar  (-). 

The  same  procedure  as  before  leads  to  the  "shot  fac- 
tors" g.  in  Table  XII  and  the  attenuation  curve  is  shown  in  Fig. 15. 

J 
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The  slope  A of  the  straight  line  and  the  mean  deviation  m of  the 
data  points  from  this  line  are 

A - -3.43792  , m=  26.7 % (8.4) 

The  scatter  is  about  the  same  as  in  the  case  of  the  first  data, 
the  slope  A is  also  about  the  same,  but  in  both  cases  are  the 
values  (8.3)  and  (8.4)  much  greater  than  what  was  obtained  in  (8.2) 
and  (8.1).  All  cases  give  greater  values  for  A than  what  is  exhi- 
bited in  Fig. 13. 


TABLE  XIII 


The  ground 

shock  factor  F 

u 

Event 

From  [2] 

From  XI 

From  [2] 

From  XII 

From  [2] 

...  . . _j 

From  XI 

1 

. 100 

.156 

. 127 

.292 

. 142 

. 322 

2 

.247 

. 128 

.315 

.351 

.264 

3 

.398 

.186 

.507 

.339 

. 566 

.384 

4 

.566 

.298 

.722 

.600 

.805 

.615 

5 

.703 

.486 

.896 

1.000 

1.000 

1.000 

6 

.785 

1.000 

1.000 

7 

1.000 

Finally  the  ground  shock  factor  F is  found  from  Fig. 13 
and  given  in  Table  XIII  with  Event  7 as  basis.  If  Event  6 is  used 
as  basis,  the  factors  will  change  as  shown  and  compared  with  the 

factors  which  can  be  computed  from  the  8 .-values  of  Table  XI.  With 

J 

Event  5 as  basis,  it  is  possible  to  compare  the  factors  from  [2] 

with  those  computed  from  the  B .-values  of  both  Table  XI  and  XII. 

J 

Even  though  the  latter  two  agree  reasonably  well,  it  must  be  con- 
cluded that  the  data  upon  which  the  present  investigation  is  based 
do  not  give  results  which  conform  with  those  reported  earlier  [2] 
from  the  same  experiment. 

9.  Final  remarks 


The  CENSE  experiment  distinguishes  itself  from  the 
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VIKERSUND  experiment  in  the  basic  philosophy  upon  which  it  is 
planned.  The  CENSE  experiment  considers  cases  where  the  detonating 
charge  may  be  placed  al  all  possible  levels  above  as  well  as  below 
the  ground  surface.  It  is  assumed  that  the  shock  waves  created  in 
these  cases  can  be  related  to  each  other  and/or  to  the  fully  con- 
tained case  by  a ground  shock  factor. 

The  VIKERSUND  experiment  is  based  on  the  philosophy 
that  a physically  important  difference  exists  between  cases  in 
which  the  shock  wave  is  created  in  connection  with  a crater  formation 
and  cases  where  this  is  not  the  case.  The  VIKERSUND  experiment 
thus  concentrates  entirely  on  cases  where  the  charge  is  detona- 
ted above  the  surface  or  in  contact  with  it.  The  case  of  the  socalled 
half-burried  charge  is  considered  separately  in  [3],  and  in  spite 
of  the  inadequate  experimental  support,  it  is  believed  that  a way 
is  suggested  whereby  the  properties  of  the  rock  both  as  a wave- 
transmitting  medium  and  its  ability  to  withstand  cratering  is  taken 
into  account. 


Because  of  this  situation,  the  two  experiments  con- 
sidered in  this  report  do  not  cover  the  same  physical  situations. 
It  could  beforehand  be  expected,  that  they  might  in  some  way  com- 
pliment each  other,  but  the  discrepancies  discovered  prevent  that. 

The  final  result  as  far  as  the  VIKERSUND  experiment 
is  concerned  is  thus  contained  in  the  two  statements  (7)  and 
equations  (6.1)  and  (6.3).  It  should  be  emphasized,  that  further 
experimentation  seems  necessary  to  furnish  an  adequate  experimen- 
tal support  for  these  results. 
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